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The energetics of mixing and defect ordering in solid solutions of fluorite-structured Th02 with 
oxides of trivalent cations (Sc, In, Y, Nd, La) are investigated by electronic density-functional- 
theory (DFT). Through DFT calculations of structures enumerated by lattice-algebra techniques, 
we identify the lowest-energy patterns for defect clustering for four separate dopant concentrations. 
The most stable structures are characterized by a repulsive interaction between nearest-neighbor 
vacancies on the oxygen sublattice. The enthalpies of formation with respect to constituent oxides 
are positive for all dopants considered, and show a tendency to decrease in magnitude as the size and 
electronegativity of the trivalent dopant decrease. Due to the small positive formation enthalpies and 
low oxygen-vacancy binding energy with La dopants, La203-Th02 solid solutions are predicted to 
have relatively high ionic conductivities relative to those for the other aliovalent dopants considered. 
Our results are compared with those for the more widely studied Zr02 and CeC>2 fluorite-structured 
solid solutions with trivalent cations. 



I. INTRODUCTION 

Design of solid-state oxygen ion conductors, oxygen 
sensors and separation membranes^^ typically involves 
aliovalent cation doping of metal oxides to form charge- 
compensating oxygen vacancies. The observation of 
high oxygen conductivity in aliovalently doped fluorite- 
structured oxides has motivated extensive experimen- 
tal and theoretical research aimed at relating the struc- 
tural and chemical defect properties in these systems to 
their ionic conductivities. Recent studies of trivalent- 
doped zirconia, hafnia, ceria and thorigPIIIl have revealed 
a strong correlation between oxygen conductivity and 
the mixing thermodynamic properties. Specifically, en- 
thalpies of formation (AH*) derived from calorimetric 
measurements tend to become more positive with dopant 
concentration (x) at low solute concentrations, in a man- 
ner consistent with a regular solution behavior; near 
some maximum composition (x max ), however, AH? is 
observ ed to le vel off or become less positive with further 
dopingP 11 | 12 | The value of x max is found to coincide with 
the maximum in measured ionic conductivity.^ These 



experimental results suggest that the mobile oxygen va- 
cancies begin to strongly associate and form stable, rela- 
tivley immobile, cl usters f or dopant concentrations near 
and beyond x rooa ,P® 12|ia | However, the detailed nature 
of these defect clusters, and how they vary in different 
systems, remains incompletely understood. 

In the technologically important yttria-stabilized zirco- 
nia (YSZ) system, atomistic simulations ^^ a gree with 
neutron and X-ray diffraction experiments^^ that oxy- 
gen vacancies tend to align along the <111> direction as 
third nearest neighbors on the oxygen sublattice, and pre- 
fer to associate with the smaller host Zr ions rath er than 
with the larger Y dopants. Further, DFT studies^HSUHl 
have shown that these defect-clustering tendencies result 
in the stability of two ordered ground-state structures: 
the experimentally observed <5-Zr3Y40i2 compound and 
an as-yet unobserved Zr02-rich phase. Through studies 
of other trivalent and divalent dopants in ZrC>2, it has 
been shown that the formation enthalpies of doped Zr02 
compounds are most positive for small dopant ions. The 
results also demonstrate a tendency for oxygen vacancies 
to associate with the smaller cation species (Zr or dopant) 



2 



in a given system. In classical atomistic simulations of 
the energetics of defec t cluster formation in trivalent - 
doped Ce0 2 (Ref! 21 * 22 1 ) it has been demonstrated that 
the atomic geometry of defect clusters is controlled by 
a combination of electrostatic and elastic energy effects, 
the latter being dependent on the dopant ion size. The 
binding energies between dopant cation and oxygen va- 
cancies also have been shown to be a strong function of 
the dopant ion size, with the smallest dopants exhibiting 
greatest binding energies with nearest-neighbor oxygen 
vacancies. 

Solid solutions of TI1O2 with aliovalent trivalent 
cations also can be utilized for solid electrolyte appli- 
cations, similar to doped zirconia and ceria. ThC>2 also 
has been gaining increasing attention due to the renewed 
interest in the Th-U fuel cycle which offers a number 
of advantages over the U-Pu cycle including: higher 
thermal conductivity of ThC>2 relative to UO2 based fu- 
els and lower concentrations of long-lived transuranic 
elements! 23 * 24 ! The trivalent dopants Y, Nd and La are 
important products in the thorium-based nuclear cycle. 
Recently, Y 2 3 - and La 2 3 -doped Th0 2 have been the 
subject of calorimetric studies of mixing energetics.^ As 
in doped CeC>2 systems, the oxygen-ion conductivity in 
Y2C>3-doped ThC>2 shows the aforementioned correlation 
between the maximum in formation enthalpy and oxygen 
conductivity. 

The key to understanding the thermodynamics of de- 
fect clustering and its correlation with conductivity is 
to identify underlying interatomic interactions between 
dopants and oxygen vacancies. In solid solutions of alio- 
valently doped oxides, electrostatic interactions between 
ions with different charges and strain-mediated interac- 
tions due to the size mismatch between host and dopant 
atoms, and the presence of oxygen vacancies, are general 
in nature. Thus, if these are the primary factors driving 
defect association, the nature of the defect clusters might 
be expected to be generic across fluorite-structured sys- 
tems, depending only on the charge and size of the aliova- 
lent dopant relative to the host. The analysis of driving 
forces for defect clustering can be more complex if other 
factors play a significant role . For e xample, in the case 
of YSZ, it has been suggeste d 14 * 25 * 26 ! that the preference 
of Zr cations to be seven-fold coordinated by oxygen in 
the lowest-energy monoclinic phase of Zr0 2 leads to a 
driving force for association with oxygen vacancies in the 
eight-fold coordinated fluorite structure. 

TI1O2 is an interesting system from the standpoint of 
studying defect interactions in fluorite-structured com- 
pounds, as (in contrast to Zr02) it is a system that dis- 
plays a stable fluorite structure in pure form, and (in 
contrast to Ce02) it displays a low equilibrium solubil- 
ity of oxygen vacancies and no tendency toward reduc- 
tion. In this study, we exploit first-principles DFT-based 
calculations combined with the lattice-algebra structure 
enumeration technique to study mixtures of ThC>2 with 
trivalent dopants. To investigate the energetics and the 
nature of defect clustering, we examine a set of dopants 



TABLE I. The effective ionic radii of ShannorP r (A) for 
trivalent dopants in the octahedral coordinatiorQand Pauling 
electronegativities X- 





Sc 


In 


Y 


Nd 


La 


r 


0.75 


0.80 


0.90 


0.98 


1.03 


X 


1.36 


1.78 


1.22 


1.14 


1.10 



a The ionic radius of the host Th 4+ is 0.94 A and its Pauling 
electronegativity is 1.3. 



spanning a range of ionic sizes and electronegativities, as 
shown in Table [H 



II. COMPUTATIONAL METHODOLOGY 

First-principles calculations were carried out within 
the DFT plane wave formalism using the GGA-PBE 
exchange-correlation functional and the projector aug- 
mented wave (PAW) potentials^ as implemente d in t he 
Vienna Ab initio Simulation Package (VASP)P^ ft 
has been demonstrated 3 ^ that the inclusion of Hubbard- 
type on-site electron interaction within the DFT + 
U approach 3 ^ for calculation of structural, elastic and 
electronic properties of TI1O2 changes the results only 
marginally as compared with pure DFT calculations. We 
employ the potentials labeled Th, O, Sc_sv, In_d, Y_sv, 
Nd and La in the VASP PBE library; for Sc_sv, Y_sv 
and In_d these potentials include semicore states as va- 
lence electrons. In all calculations we employ a cutoff 
energy of 500 eV which is 25% higher than the default 
value for the potentials chosen to ensure convergence 
of the calculated stress tensors for volume relaxations. 
All atomic structure relaxations were performed start- 
ing from configurations of atoms in ideal fluorite posi- 
tions and optimizing cell volume, cell shape and all ionic 
positions within a cell by applying a conjugate-gradient 
algorithm until atomic forces were converged to a mag- 
nitude less than 0.1 meV/A. In these relaxation runs, we 
employed the Monkhorst-Pack scheme^ to sample the 
Brillouin zone, using meshes with a density of at least 
4x4x4. Subsequent to the structural optimizations, an 
additional static run for each final optimized geometry 
was made using accurate tetrahedron integration with 
Blochl corrections,^ with a sufficiently dense sampling of 
the Brillouin zone to converge the total energy to better 
than 10 -6 eV. For structures involving the trivalent Nd 
dopant, which contains three valence / electrons, calcu- 
lations were performed spin-polarized, using a ferromag- 
netic arrangement of the local magnetic moments (for a 
few structures both ferromagnetic and antiferromagnetic 
ordering was considered, and the former was found to be 
lower in energy). 

In fluorite-structured ThC>2 doped by oxides of triva- 
lent cations, one has to account for binary disorder on 
both cation (Th 4+ , M 3+ , where M denotes the trivalent 
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FIG. 1. TI1O2 cubic fluorite structure with two triva- 
lent dopants on the cation fee sublattice and one charge- 
compensating oxygen vacancy on an anion simple cubic sub- 
lattice. 



TABLE II. Formation enthalpies (kj/mol-cation) of the most 
stable defect associates for Thi-^M^O^-o.Sx (M = Sc, Y, 
La) systems, as derived from calculations employing 2x2x2 
supercells of the host TI1O2 conventional unit cell. The stable 
defect clusters correspond to a nearest-neighbor M-0„ ac -M 
triangular cluster. 



Sc 



Y 



0.125 
0.0625 



0.0625 



-1.9 (-1.9) 



-1.1 (-1.1) 



La 



Formation enthalpy 

14.4 6.3 3.4 

8.8 3.4 2.0 

Binding energy 



-0.2 (-0.6) 



dopant species) and anion (0 2 ~, O vac ) sublattices (see 
Fig. [T]). In order to probe the configurational space and 
sample different possible atomic arrangements for given 
dopant concentrations, we employed a structure enumer- 
ation technique making use of the multicomponent multi- 
sublattice formalism previo usly dev eloped to study alloys 
and recently applied to YSZpEEU We performed a com- 
plete structure enumeration in supercells consisting of up 
to nine primitive fee unit cells (for the present study we 
focus on compositions where we enumerated structures 
only up to six unit cell s) usi ng the Alloy Theoretic Auto- 
matic Toolkit (ATAT)PHS These enumerated structures 
were considered for dopant concentrations of a;=0.33 and 
0.5 cation mole fraction. We also employed supercells, 
composed of 2x2x2 conventional fluorite unit cells (96 
atoms), to study defect energies at the more dilute com- 
positions of £=0.0625 and 0.125, corresponding respec- 
tively to a single vacancy with two dopant cations, and 
two vacancies with four dopants. 

Formation enthalpies of solid solutions normalized per 
cation can be expressed with respect to the constituent 
oxides as 



AH* = E[Th 1 - x M x 2 -o.6x]-(l-x)E[Th0 2 ]-xE[M0 1 . B ], 



where E\Th0 2 \ denotes the energy (per cation) of 
fluorite-structured TI1O2. Similarly, E[MOt,s] denotes 
the energy (per cation) of SC2O3, In 2 03 and Y2O3 in the 
experimentally-observed cubic bixbyite (C-type sesquiox- 
ide) structure, or Nd2C>3 and La2C>3 in the observed low- 
temperature hexagonal A-type structure. The energy per 
cation for a solid solution with a cation dopant concen- 
tration x is denoted as E[Thi- x M x 2 _ x / 2 ]. With the 
above definition, positive values of AH* correspond to 
an energetic tendency for phase separation of the solid 
solutions into the respective constituent oxides. 



III. RESULTS 

TI1O2 compounds form in the cubic fluorite structure 
in which the cations form an fee sublattice and anions 
form a simple cubic sublattice (Fig. [IJ. Doping TI1O2 
by trivalent cations introduces one charge-compensating 
oxygen vacancy per two trivalent dopants and leads to 
a coupled intra- and intersublattice interactions between 
all the species present (Th 4+ , M 3+ , O 2 ", O vac ). To in- 
vestigate the nature of these interactions, and the ge- 
ometries of the defect clusters that they favor energeti- 
cally, we investigate the formation energetics of a variety 
of superstructures for relatively dilute and concentrated 
compositions in sub-sections A and B, respectively. 

A. Defect association in the dilute regime 

In this section, we focus on systems with relatively di- 
lute dopant concentrations employing a 2x2x2 supercell 
of the host TI1O2 crystal with one and two electrically 
neutral M-0„ ac -M units, corresponding to x=6.25% and 
12.5% dopant concentrations, respectively. We focus 
on a comparison of the energetics for under-sized Sc, 
intermediately-sized Y and over-sized La dopants (Ta- 
ble I). We examine a number of plausible defect arrange- 
ments and obtain only positive formation enthalpies in 
agreement with recent calorimetric data for Y and La 
dopants 

For the most dilute dopant concentration (6.25%) the 
most energetically stable clusters correspond to the com- 
pact triangle geometry involving M and O vac as nearest- 
neighbor pairs. The calculated formation energies of 
these compact clusters are listed in Table [TTJ and are 
seen to significantly decrease in magnitude from Sc to 
Y to La. We also study the energetics of clustering be- 
tween two triplet defect clusters (12.5% defect level) by 
considering a number of different atomic arrangements. 
We find that the lowest-energy arrangements correspond 
to two nearest-neighbor triplet M-0„ ac -M clusters, with 
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the two oxygen vacancies aligned as second neighbors on 
the oxygen sublattice, i.e., along the <110> direction, 
with three dopant cations and one host thorium ion com- 
posing tetrahedron around the oxygen vacancy. The for- 
mation energies of these clusters are also listed in Table 
[TT| and again show a trend of decreasing magnitudes for 
AH? going from Sc to Y to La. 



We further estimate the binding energies between the 
trivalent dopants and oxygen vacancies for the most sta- 
ble defect clusters. The binding energy here is defined as 
the difference in energy between a system with a single 
defect cluster and one with isolated defects, such that 
negative values indicate an energetically stable defect 
cluster. Due to the relatively small supercells used in 
this study, we cannot derive highly accurate values for 
these binding energies. Rather we explore trends by es- 
timating the energy of the isolated defects as the energy 
of the supercell in which the defects are separated as far 
as possible within the 2x2x2 supercell. 



In the case of the M-Oi, ac -M triplet cluster, the 
binding energy can be split into two parts correspond- 
ing to subsequent breaking of the first and the second 
nearest-neighbor M-0„ as bond. This allows one to es- 
timate whether or not the dopant-vacancy interactions 
are pairwise-additive in nature. The energy to break the 
first and second nearest-neighbor M-0„ ac bonds is given 
in the final row of Table [TTJ For Sc and Y the interac- 
tions are found to be relatively short-range in nature, 
and the binding energies presented in Table |n] do in- 
dicate nearly pairwise-additive interactions between the 
dopant and oxygen vacancy (i.e., the energy to break the 
first and second nearest-neighbor yi-O vac bond is iden- 
tical within the precision of the DFT calculations). In 
the case of La, however, the interactions are found to 
be much longer ranged; the interaction energy between 
a La ion and a La-0„ ac pair is found to vary apprecia- 
bly with distance out to the largest separation attainable 
in the 2x2x2 supercell. For the La dopant, it can be 
seen that the interactions between dopant and oxygen 
vacancy are not pairwise additive, as it requires nearly 
three times the energy to break the second La-0„ ac bond 
as it does the first. The longer-range an non-pairwise na- 
ture of the interactions between La and O vac is likely a 
manifestation of elastically-mediated interactions, as su- 
percell calculations performed by constraining the ions 
to ideal fluorite lattice sites yielded interaction ranges 
and pairwise-additive behavior very similar to that de- 
scribed above for the Y dopants. The fact that elastic 
interactions can lead to long-range and non-pairwise in- 
teractions has been demonstrated by Bugaev et alpS It 
is nevertheless important to point out that the binding 
energies are particularly small for La system. More gen- 
erally, we note that the results in Table [TT] show a trend 
towards decreasing magnitude of the binding energies go- 
ing from Sc to Y to La. 
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FIG. 2. Formation enthalpies (kj/mol-cation) of TI14M2O11 
(M = Sc, Y and La) solid solutions corresponding to £=0.33 
doping level. Red boxes represent results obtained considering 
only volume relaxations, with the ionic positions for the end- 
member compounds and mixtures fixed at ideal fluorite lattice 
sites. The blue circles represent results obtained for fully- 
relaxed structures, and with end-member M2O3 compounds 
in their experimentally-observed low-temperature structure. 
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FIG. 3. The same as in Fig. [2]for Th 2 M 2 7 (M = Sc, Y and 
La) solid solutions corresponding to £=0.5 doping level. 



B. Defect association in the concentrated regime 

In this section we analyze the energetics of compounds 
with x=33% and 50% dopant concentrations. Both com- 
positions are beyond the dopant levels x max where the 
maximum in the measured formation enthalpies are ob- 
served for Y and La dopants; 3 ^ these concentrations are 
thus expected to be in the regime where the dopant 
cations and oxygen vacancies interact strongly. The com- 
position A2B2O7 (50% dopant level) also is associated 
with the widespread group of pyrochlore-structured crys- 
tals. The pyrochlore structure is a derivative of fluorite, 
in which A and B cations are ordered along the <101> 
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direction, and oxygen vacancies reside in the tetrahedral 
interstices between adjacent B cations. 

Compared to the dilute concentrations discussed 
above, the number of distinct ionic configurations are 
much larger at the higher concentrations. In order to 
determine the energetically preferred atomic arrange- 
ments, we performed an exhaustive structure enumer- 
ation for Th 4 M 2 0ii (33% dopant level) and Th 2 M 2 7 
(50% dopant level) compositions. The former enumera- 
tion yielded 117 symmetry-distinct configurations with 
unit cells containing 4 Th and 2 M cations, 11 oxy- 
gen ions and 1 oxygen vacancy, arranged over the sites 
of the ideal fluorite structure. Similarly, the enumera- 
tion at 50% dopant concentration yielded 27 symmetry- 
distinct fluorite superstructures containing 2 Th and 2 
M cations, 7 O ions and 1 oxygen vacancy per unit cell. 
The calculated energies for each of these structures are 
plotted in Figs. [2] and [3] for £=0.33 and 0.5, respec- 
tively. In these plots, results labeled "relaxed" corre- 
spond to formation energies computed as the energy dif- 
ference between a fully relaxed Thi_ 2; M 2 ,02-o.5x super- 
structure and the concentration-weighted average of the 
energies of relaxed fluorite-structured Th0 2 and M 2 C>3 
in its experimentally-observed low-temperature crystal 
structure. The "unrelaxed" results correspond to the 
analogous energy difference, but using structures that are 
only relaxed with respect to volume, constraining all of 
the ions to occupy the sites of an ideal fluorite lattice. 

The results in Figs. [2] and [3] show a large spread in 
the DFT-calculated formation enthalpies at each fixed 
composition. For the relaxed structures this spread in 
AH? between the lowest and highest energy structures 
is largest for Sc (85.5 and 87.8 kJ/mol-cation at £=0.33 
and 0.5, respectively) and smallest for La (56.6 and 39.6 
kJ/mol-cation at x=0.33 and 0.5, respectively). All of 
the formation energies for the relaxed structures are pos- 
itive in Figs. [2] and [3j indicating a preference for phase 
separation at low temperatures, consistent with experi- 
mental measurements for La and YP^J Further, the cal- 
culated values of AH? decrease in magnitude going from 
Sc to Y to La. This trend is consistent with the results 
discussed above for the more dilute compositions, and 
with measured solubility limits at T=1773 K, which are 
around 50 mol% for La and 20 mol% for Yp^ and negli- 
gible in the case of Sc (the solubility of Sc has not been 
measured to the best of our knowledge, but is expected 
to be much smaller than for La and Y). 

The differences in AH? between the relaxed and unre- 
laxed structures plotted in Figs. [2] and [3] reflect a balance 
between two competing effects: (i) the stabilization of 
the Thi_ x M a; 2 compounds by atomic relaxation, and 
(ii) the stabilization of the end- member M 2 03 compound 
when it is allowed to change from an ideal fluorite lattice 
to its experimentally observed low-temperature crystal 
structure. The first and second effects give rise to con- 
tributions that decrease and increase the value of AH? , 
respectively. For the Y dopant, with a size that is most 
closely matched to Th, the first and second effects are 




FIG. 4. Illustration of the stable defect ordering pattern in 
Thi-zM^O^-o.Sa; solid solutions. The structure is character- 
ized by charge-neutral tetrahedra composed of one trivalent 
dopant M, two corner-shared M and one Th ion with the oxy- 
gen vacancies aligned as second nearest neighbors along the 
<110> direction. 



comparable in magnitude for the low-energy structures 
at each composition, such that the formation energies are 
nearly identical for relaxed and unrelaxed structures; the 
first effect is considerably larger than the second for the 
high-energy structures, leading to a decrease in the mag- 
nitude of AH^ for the relaxed versus unrelaxed struc- 
tures. For the under-sized Sc dopant, relaxation has the 
effect of lowering the formation energy for all structures 
considered, while for the over-sized La ion the low-energy 
and high-energy structures are respectively destabilized 
and stabilized by atomic relaxations. It is interesting 
to note that in a recent analysis of measured formation 
energies for La-doped Th0 2 solid solutions^' the mix- 
ing energy, defined as the formation energy of the solid- 
solution relative to fluorite-structured Th0 2 and fluorite- 
structured La 2 03, was computed to be negative, even 
though the measured values of AH ' relative to fluorite- 
structure Th0 2 and A- type La 2 3 , are positiveP^J These 
results are qualitatively consistent with the present com- 
puted energies for La-doped Th0 2 , which show a change 
in the sign of the formation energy for relaxed versus 
unrelaxed low-energy structures, driven by the energy 
difference between fluorite and A-type La 2 C>3 structures. 

We consider next the qualitative features of the low 
and high-energy ionic configurations enumerated in this 
study. Table |III| summarizes the formation enthalpies of 
the most stable Th4M 2 0n and Th 2 M 2 07 compounds, 
along with their main structural features, including the 
direction connecting oxygen-vacancy neighbors, and the 
number of cations of each type (Th and M) in the nearest- 
neighbor tetrahedron surrounding the vacancy. We in- 
clude such information for several of the lowest-energy 
structures for each dopant; some of the structures have 
nearly degenerate or close energies and therefore a num- 
ber of defect ordering patterns may be observed particu- 
larly in the presence of entropic effects at elevated tem- 
peratures. 

We observe that the oxygen vacancies prefer to align as 
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TABLE III. Formation enthalpies of Thi-^M^Cb-o.Sx (M = 
Sc, Y, La) structures are listed for the lowest-energy fully- 
relaxed structures enumerated in this study. The structures 
listed have formation energies within 2 kj/mol-cation of the 
lowest energy structure for each dopant at each composition 
considered. The third column lists the direction of ordering 
of oxygen- vacancy chains. The last column lists the num- 
ber of cations of host (Th) and dopant (M) species in the 4 
tetrahedral nearest-neighbor positions surrounding an oxygen 
vacancy. 

Nearest Neighbors 



.X' 


M 


AH f 


Direction 


Th, M 


0.5 


Sc 


28.1 


^ A- -LKJ ^ 


1 3 






28.3 




3 1 




Y 


14.8 


<110> 


1 3 




La. 


9.1 


<110> 


1, 3 






9.2 


<110> 


2, 2 






10.1 


<110> 


1, 3 






10.5 


<110> 


2, 2 


0.33 


Sc 


22.9 


<110> 


1, 3 






22.9 


<110> 


3, 1 






23.4 


<110> 


1, 3 




Y 


11.7 


<110> 


1, 3 






13.0 


<110> 


1, 3 






13.2 


<110> 


1, 3 






13.5 


<110> 


2, 2 






13.6 


<110> 


1, 3 




La 


4.4 


<112> 


4, 






6.1 


<112> 


2, 2 






6.2 


<112> 


3, 1 



second nearest neighbors along the <110> directions for 
all dopants with an exception of La for a;=0.33 where the 
vacancies are ordered along the <112> direction. The 
ordering of oxygen vacancies along the <110> direction is 
also observed for the La system at x=0.33, but at slightly 
higher energies (with AH f > 6.7 kJ/mol-cation). An 
analysis of nearest neighbors around the oxygen vacancies 
shows that the only case in which the oxygen vacancy is 
surrounded completely by four host Th cations is for the 
lowest-energy La-doped TI1O2 structure at x=0.33. In all 
other cases we see that the tetrahedra of cations around 
the vacancy, as pictorially shown in Fig. [4j are composed 
of both host Th and dopant M ions, with a slight overall 
preference for binding to trivalent dopant species. 

For Y dopants at x=0.33 and 0.5, the lowest-energy 
structures are shown in a <110> projection in Fig. [7] In 
this projection, the smaller (red) circles denote <110> 
rows of oxygen ions containing no vacancies, while the 
darker and lighter larger (blue) circles denote <110> 
columns of Y and Th cations, respectively; the black open 
boxes in the oxygen rows correspond to <110> rows of 
oxygen vacancies. Both structures in Fig. [7] contain the 
same structural motif, which, in the projected image in- 
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volves a column of oxygen vacancies located between a 
triangle of two Y and one Th ion. This motif is illustrated 
in a different representation in Fig. [4] where it is more 
clearly apparent that the lowest-energy Y-doped ThC>2 
structures are characterized by a clustering of vacancies 
as second neighbors on the anion sublattice (i.e., along 
<110>) directions, and a clear association of the oxy- 
gen vacancies with nearest-neighbor dopant cations. In 
Fig. [4] the motif of the defect cluster involves tetrahedra 
consisting of three Y and one Th cation, surrounding an 
oxygen vacancy; since these tetrahedra are arranged in 
a corner-sharing configuration the defect clusters are net 
charge neutral. For the TI12Y2O7 structure (x=0.5) the 
lowest-energy structure shown in Fig. 7] (left) can also be 
identified as being characteristic of the pyroclore struc- 
ture with the oxygen vacancies residing in the tetrahedral 
interstice between adjacent Y dopants, aligned along the 
<110> direction. 

The structures of the highest-energy configurations 
for Y-doped ThC>2 considered in our calculations are 
illustrated in Fig. [HJ These structures are character- 
ized by the location of rows of oxygen vacancies be- 
tween host Th ions, and thus contain large regions that 
are locally charge imbalanced. The large energy differ- 
ence between these configurations and those characteris- 
tic of the lowest-energy structures demonstrates a strong 
preference for association of oxygen-vacancy chains with 
dopant cation species in this system. 



IV. DISCUSSION 



The results in Figs. [2] and [3] and Tables [TT] and |TTT] fea- 
ture a strong variation in the calculated values of AH? as 
the dopant species change from Sc to Y to La. To gain 
further insight to the origins of this trend, we explore 
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the correlations between AH? and the dopant properties 
listed in Table [Tj Figure [5] plots calculated formation en- 
ergies versus ionic radius for the lowest energy structures 
containing Sc, Y, and La dopants, as well as the two addi- 
tional trivalent dopants In and Nd. A tendency towards 
decreasing magnitudes of formation energies for dopants 
with increasing size is observed, but with the notable ex- 
ception of In which has the largest positive calculated 
value of AH? and an ionic radius larger than Sc. 

An interesting correlation is shown in Fig. [6j which 
plots AH$ versus the Pauling electronegativity (%) of 
the dopant species. A monotonic decrease is observed in 
the calculated values of AH* with decreasing %. Since 
a decrease in the electronegativity of the dopant cation 
corresponds to an increasing electronegativity difference 
between it and oxygen, the trend in Fig.[6]suggests that a 
decreasing ionic character in the bonding between dopant 
and oxygen ions has a tendency to destabilize doped 
Th02 mixtures. 

To investigate this correlation further, we performed 
an analysis of the electronic structures of the lowest- 
energy enumerated compounds, through calculations of 
Bader charges and partial electronic densities of states 
(DOS). For low-energy Sc and La doped structures, the 
calculated Bader charge on La is larger than that on Sc 
by appoximately 0.1 electrons, consistent with the higher 
value of x f°r the latter. The change in the ionic char- 
acter of the bonding between oxygen and dopant cations 
with varying \ is a l so indicated in the calculated par- 
tial DOS. The calculated DOS feature a partial covalent 
character to the bonding, characterized by a weak hy- 
bridization of O 2p and dopant d electrons. The degree 
of hybridization, as measured from the integrated num- 
ber of electrons in the O 2p band, is found to decrease 
monotonically going from Sc to In to Y to Nd to La, con- 
sistent with the trend in the dopant electronegativities. 



Overall, the analysis suggests that the mixing energetics 
of aliovalently-doped TI1O2 solid-solutions is influenced 
not only by the size and charge of the dopant species, 
but also by chemical effects associated with the degree 
of covalency characterizing the bonding between dopant 
and oxygen ions. 



It is interesting to compare the results of the present 
study with those obtained in previous investigations of re- 
lated trivalent-doped fluorite-structured solid solutions. 
In the most stable <5-phase of YSZ each oxygen vacancy is 
enclosed in a tetrahedron consisting of two Y and two Zr 
cations, and the vacancies are ordered along the <111> 
direction, as third-nearest neighbors. This atomic ar- 
rangement exhibits a negative formation enthalpy, in- 
dicating a compound-forming tendency. On the con- 
trary, in trivalently-doped TI1O2 we observe only posi- 
tive enthalpies of formation for all dopants considered, 
and an overall preference for <110> ordering of oxy- 
gen vacancies (i.e., as second-nearest neighbors on the 
anion simple-cubic sublattice). In agreement with the 
present results, it has been observed that formation en- 
thalpies of aliovalently-doped Zr02 are most positive for 
small dopants The same trend has been found for 
pyrochlore-structured solid solutions of Ce02 with triva- 
lent dopants pi where the formation energy increases in 
magnitude going from La to Gd to Y dopants. It also 
has been demonstrated in the dilute regime that small 
dopants exhibit the strongest binding with the neighbor- 
ing oxygen vacancies in Ce02, which is again consistent 
with the results of the present study. 



The variations in binding and formation energies found 
for the different dopants considered in this study may 
have important consequences for oxygen ionic conductiv- 
ity in aliovalently-doped TI1O2 compounds. Specifically, 
La dopants are found to have significantly smaller forma- 
tion energies and binding energies (see Tabic [Tl| than the 
other dopants considered in this work. Further, we found 
that the La-doped compounds showed a higher degener- 
acy in the energy of competing structures at more con- 
centrated compositions. Overall, these trends are con- 
sistent with experimental observations. Specifically, the 
solubility of La in Th02 is significantly larger that that 
for Y or Sc doped systems, and from an analysis of the 
calorimetry results in RefP^the association energy for de- 
fects in La-doped TI1O2 is estimated to be considerably 
smaller than in Y-doped TI1O2. The smaller magnitudes 
for the formatin energies and oxygen- vacancy binding en- 
ergies for La dopants is expected to give rise to higher 
ionic conductivities in this system, relative to Y-doped 
TI1O2. While conductivity data for the latter system are 
available,^ measurements for the former have not yet 
been published, to the best of our knowledge. Such mea- 
surements for La-doped Th02 would thus be interesting 
to pursue in light of the present results. 
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FIG. 7. Lowest-energy atomic arrangements in Thi-^Y^C^-o.Sz: solid solutions at £=0.5 (left) and 0.33 (right). The different 
circle symbols correspond to rows of O, Th and Y ions directed along the <110> direction (out of the plane of the page), while 
the open boxes correspond to rows of oxygen vacancies along the same direction. For clarity all atoms are displayed in their 
ideal fluorite positions. 
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FIG. 8. Highest-energy atomic arrangements in 

Thi_ :E M a ;02-o.5i solid solutions at x=0.5 (left) and 
0.33 (right). The symbols have the same meanining as in 
Fig. [7] For clarity all atoms are displayed in their perfect 
fluorite positions. 



V. CONCLUSIONS 

In this study we have employed DFT calculations to 
study formation energetics and defect-ordering tenden- 
cies in Th02 compounds doped with the trivalent cation 
species: Sc, In, Y, Nd and La. The main results derived 
from these calculations can be summarized as follows: 

(1) Solid solutions of fluorite-structured ThC>2 with ox- 
ides of trivalent cations (Sc, In, Y, Nd, La) exhibit pos- 
itive enthalpies of formation from constituent oxides, in- 
dicative of a phase-separating tendency across all dopant 
concentrations considered (6.25, 12.5, 33.3 and 50%). 
This strongly suggests that there are no underlying stable 
long-range ordered phases in these compounds. However, 
the large variation in energy for different dopant / vacancy 
arrangements suggests an appreciable driving force for 
short-ranged ordering (clustering) in these solid solutions 
at finite temperatures. 

(2) For a given dopant concentration the formation 
enthalpies tend to decrease in magnitude as the size and 
electronegativity of the trivalent dopant decrease. This 



trend is consistent with available experimental data on 
solubility limits which are highest for Nd- and La-doped 
Th0 2 . 

(3) The lowest-energy atomic structures are character- 
ized by a repulsion of oxygen vacancies as nearest neigh- 
bors. Oxygen vacancies are found to prefer to align as 
second nearest neighbors along the <110> direction for 
all of the dopants and concentrations considered, with 
the exception of 33.3% La-doped Th02 where they order 
as more distant neighbors along the <112> direction. 

(4) The solute-vacancy binding energies in Th0 2 solid 
solutions with under-sized Sc and intermediate-sized Y 
dopants are estimated to be on the order of a few kJ /mol- 
cation. For over-sized La dopants, this binding energy is 
estimated to be significantly smaller in magnitude. 

(5) The small magnitude of the formation and binding 
enthalpies for La-doped Th0 2 solid solutions is expected 
to result in an ionic conductivity that is higher in this 
system than for the other trivalent-doped systems con- 
sidered in this study. 
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